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Stability of contact discontinuity for the 
Navier-Stokes-Poisson system with free boundary 

Shuangqian Liu* Haiyan Yin} Changjiang Zhu* 


Abstract 

This paper is concerned with the study of the nonlinear stability of the contact discontinuity 
of the Navier-Stokes-Poisson system with free boundary in the case where the electron background 
density satisfies an analogue of the Boltzmann relation. We especially allow that the electric potential 
can take distinct constant states at boundary. On account of the quasineutral assumption, we first 
construct a viscous contact wave through the quasineutral Euler equations, and then prove that such 
a non-trivial profile is time-asymptotically stable under small perturbations for the corresponding 
initial boundary value problem of the Navier-Stokes-Poisson system. The analysis is based on the 
techniques developed in HD and an elementary L 2 energy method. 
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1 Introduction 

1.1 The problem 

The dynamics of the charged particles in the collisional dusty plasma can be described by the Navier- 
Stokes-Poisson (denoted as NSP in the sequel) system [15,. The one-dimensional NSP system in the 
Eulerian coordinates takes the form of 

' d t p + d x (pu) = 0, 

|% ( pu ) + d x ( pu 2 + p) = pd x (j) + pd x u, 

d t W + d x (Wu +pu) = pud x § + pd x {ud x u) + nd x 9, 

= P ~ Pe(4>)- 
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The unknown functions p, u and 6 stand for the density, velocity and absolute temperature of ions, 
respectively, while p > 0 is the viscosity coefficient and k > 0 is the heat conductivity coefficient. W 
stands for the total energy of the ions, taking the following form: 



P 

7 - 1 ’ 


where 7 > 1 is the adiabatic exponent, p is the pressure which is given by 

p = RpO = Ap 1 e ' LrrS , 


where S is the entropy and A , R are both positive constants. The self-consistent electric potential 
0 = cj>(x,t) is induced by the total charges through the Poisson equation. The density p e = p e {(/>) of 
electrons in ED depends only on the potential in the sense of an analogue of the so-called Boltzmann 
relation, cf. El [ 20 ]. Specifically, through the paper we suppose that 

(.4) p e {4>) ■ (0m 1 0 m ) ( Pm, Pm ) is a smooth function with 

p m = inf p e (0), p M = sup p e (0), 

satisfying the following two assumptions: 


(4l) p e (0) = 1 with 0 G (0m, 0 m); 

(A 2 ) pe{4>) > 0, Pe(0) < 0 for each 0 G (0 m , 0 m)- 


The assumption (Ai) just means that the electron density has been normalized to be unit when the 
potential is zero, since the electric potential in ed can be up to an arbitrary constant. The sign of the 
first derivative of the function p e ( 0 ) in the assumption (A 2 ) plays a crucial role in our analysis, it is to 
be further clarified later on, see (11.161) . etc. 

An important example satisfying ( A ) can be given as 


Pe(0) 


1 



0m — 00 , 0M 


7e 

7e - 1 


( 1 . 2 ) 


_ 0 > 

with 7 e > 1 and A e > 0 being constants. Note that p e (0) —> e~~^ and 0 m —> +00 as j e —► 1 , 
which corresponds to the classical Boltzmann relation. In fact, Ol can be formally deduced from the 
momentum equation of the isentropic Euler-Poisson system for the fluid of electrons with the adiabatic 
exponent j e under the zero-limit of electron mass, namely, d x {A e pJ e ) = —p e d x (f>. 

In this paper, we consider the system (ED in the part +00 > x > x(t), where x = x{t) is a free 
boundary with the following dynamical boundary conditions 

dx(t) 

—^j-=u{x{t),t), x(0) = 0, (p- pd x u) | x=x (t)=l>-, 0(x{t),t) = 6-, <j)(x{t),t) = 0-. (1.3) 

We also assume 0 satisfies the boundary condition at far field: 

lim 0(x,t) = 0 + . (1.4) 

x —»-+oo 


The initial data is given by 

{p,u,6)(x, 0) = {p o ,u o ,0 o ){x), lim (p 0> u 0 ,0 o )(x) = (p + , u+, 8+). (1.5) 

a:—»-+oo 

Here p+ > 0, 8± > 0, p- > 0, u + and 0± are assumed to be constant states. Also, po{x) > 0 is 
supposed, so that the ions flow has no vacuum state. In addition, we of course assume 0o(20 satisfies the 
compatibility condition and 0 satisfies the the quasincutral condition at far field, i.e. 

#0(0) = d~, p e {4>+) = P+- (1-6) 
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Our main purpose concerns the large time behavior of solutions to O'- O- o and to ex¬ 

plore this, it is more convenient to use the Lagrangian coordinates. That is, consider the coordinate 
transformation: 


x=> p(y,t)dy, 

J x(t ) 

We still denote the Lagrangian coordinates by (x, t) for simplicity of notation. Noticing that 

r 

/ p(y, t)dy —> + 00 , as x —> + 00 , 

J x(t ) 

one sees that 0 . O', o and ( 11 . 51 ) can be transformed as the problem with fixed boundary in the 
form of 


" d t v — d x u = 0 , x > 0 , t > 0 , 

q 1 n 9 x (j) (d x u 

d t u + d x p = - 1 - pd. 

v 


v 

R , « ( d x u ) 2 

o t 0 + po x u = p -h KO : 


7-1 


V 


, x > 0 , t > 0 , 

, x > 0 , t > 0 , 


dj 

V 


d x (j) 


d x ( — ) = 1 - vp e {4>), x > 0 , t > 0 , 


with boundary condition 


<?r 


0 ( 0 , t) = 0 _, Ip — p - I ( 0 , t) = p~, <p(0, t) = (/>_, lim <j)(x, t) = 0 + , t > 0 , 


x —>-+oo 


and the initial data 


(v,u,9)(x,0) = (v o ,u o ,0 o )(x), x > 0, lim (v 0 , u 0 , 0q)(x) = (u+, u+, 9 + ). 

x—t+oo 

Here v = 1/p stands for the specific volume. Moreover, 


0 o ( 0 ) = 0 - and v + = 


0eW>+) 


hold according to O- 


(1.7) 


( 1 . 8 ) 

(1.9) 


1.2 Quasineutral Euler equations and contact waves 

In order to study the large time behavior of the solution [v(x,t),u(x,t),8(x,t),(j)(x,t)] to the initial 
boundary value problem clid , m and m , we expect that [v(x,t),u(x,i),8(x,t),^)(x,t)] tends time- 
asymptotically to viscous contact wave to the Riemann problem on the quasineutral Euler system 

r d t v - d x u = 0 , 


a * d x<t> 

o t u + o x p =-, 

V 

— - d t 8 + pd x u = 0 , 
7-1 

. 1/W = Pe((j>), 


( 1 . 10 ) 


with Riemann initial data given by 

[v,u,9]{x, 0 ) = 


[u_,u_, 0 _], x < 0 , 
[v+,u+,8 + \, x>0. 


( 1 . 11 ) 


According to H5] , one sees that the Riemann problem (11.101) and (11.111) admits a contact discontinuity 
solution 


[vC°,u CD ,9 CD A CD ] (x,t) = 


[v-,u-,8-,(j)-\, x < 0 , 
[v+,u+,8 + ,(j) + \, x > 0 , 


( 1 . 12 ) 
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on the condition that 


u- = u+, p-=p(v-,0-)=p + +p‘ l ’(v + )-p tl> (v-), (1.13) 

where 

r v i 

P+=p(v+,9+), <j>± =p~ 1 {l/v±) andp 0 =p^(v) = / dg. 

J B 3 Pe(Pe (p) 

On the other hand, due to the dissipation effect of the NSP system (HD , a viscous contact wave 
[v cd , u cd , 0 cd , (f > cd ] corresponding to the contact discontinuity [u 0 - 0 , 0 C ’ I) , cf> CD ] defined as (11.121) can 

be constructed as follows. We first denote p cd = p cd (v cd ,9 cd ) = Re cd . Since the quasineutral pressure 
p cd + p‘t’ for the profile [ v cd , u cd , 9 cd , <j) cd ) is expected to be almost constant, we set 


p~=p 


cd 



e 3 p'e(p 


- 1/1 


(*)) 


dg. 


(1.14) 


Noticing that < o and p' e {-) < 0, from which and (11.141) and the implicit function theorem, we see 
that there exists a differentiable function f(9 cd ) such that 


v cd = f(6 cd ), v± = f(0 ± ), 


(1.15) 


provided that \9+ — 9-\ is suitably small. Furthermore, by a direct calculation, it follows that 

f\0 cd ) = rA R T > 0- 
P (l) cd ) 2 p^ (0 cd ) 

We now rewrite the leading part of (|1.7[) 3 (the third equation of (11.71) ) as 

With (11.151) and (11.171) in hand, we further conjecture that [v cd , u cd , 9 cd ] satisfies 

{ d t v cd - 8 x u cd = 0, v cd = f(9 cd ), 

9 cd (0, t) = 9-, 9 cd (+oo,t) = 9 + , v cd (0 1 t)=v-, v cd (+oo, t) = v+. 


(1.16) 


(1.17) 


(1.18) 


By virtue of (11.181) . we obtain a nonlinear diffusion equation as follows: 


d t 9 cd 


g{9 cd ) 


( d x e cd \ 

V/(^ cd )y ’ 


9 cd (0,t) = 9-, 9 cd {+oo 1 t) = 9+, 


(1.19) 


where g{9 cd ) = —+ p cd f{9 cd ) > 0. Applying the same argument as in [l], one sees that (11.191) admits 
a unique self similarity solution 0 cd (£), £ = / X +t - Additionally, it turns out that 9 cd is a monotone 

function, increasing if 9 + > 9- and decreasing if 9 + < 9-, and more importantly, one can show that 
there exists some positive constant <5, such that for <5 = 1 9 + — 6L | < S , 9 cd satisfies 

(1 +t) \d 2 x 9 cd \ + (l + t)% \d x 9 cd \ + \9 cd -9±\ <C6e~^ : as i->+oo, (1-20) 

where c\ is some positive constant. After 9 cd and v cd are obtained, we now define [u , </> cd ] as follows 


icd 


U cd = u + — K 


= Pe( VO. 

+oo/,( 0 cd) / d x 9 cd \ 

s(0"-) *(/(« rf )J 


f 


= U+ + 


* f(9 cd ) 


-dj) 


cd 


+oo 


(' d x e cd ) 2 


g(0 cd )f(9 cd ) x J x f(9° d ) \g 
4> cd (ti,t)=<f>-, c/) cd (+oo,t) = </>+, u cd (+oo,t) = u+. 


( 1 . 21 ) 


— {9 cd )d. 
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It should be noted that </>± = p~ 1 (l/v±), and u cd (0,t) may not equal to u+. 

In view of (11.121) . (11.181) . (11.201) and (11.211) . it is straightforward to compute that [v cd , u cd , 9 cd , (f cd ] 
satisfies 

\\[v cd -V CD ,u cd -U CD ,e cd -0CD i ^ d ^ ( j ) CD]\\ Lp ^ =o(k^) (1 +t)&, p> 1, 

which implies the viscous contact wave [v cd , u cd , 9 cd , <p cd ] (x,t) constructed in (11.181) and (11.211) approx¬ 
imates the contact discontinuity solution [v CD , u CD , 0 CD , 4> CD ] to the quasincutral Euler system (11.101) 
in L p norm, p > 1 on any finite time interval as the heat conductivity coefficients k tends to zero. 
Moreover, we see that the viscous contact wave [v cd , u cd , 9 cd , <j> cd ] (x, £) solves the Navier-Stokes-Poisson 
system O) time asymptotically, that is, 


where 


d t v cd - d x u cd = o, 


cd 


d t u cd + d x p cd = ^- + pd x 


d x u cd 

n,Cd 


R aacd , cdo cd ( d x U cd ) 2 

d t 0 C +p d x U = p -—- \-K& : 


IZi, 

d T e cd 


7-1 


it cd 


nCd 


1 Z 2 , 


cd 


dx4>' 

n,Cd 


= l-v cd p e {cj> cd )+n: 3 , 


=dt 


( d 


/» + 00 


K (d x d cd ) 2 ffY 


\g(0 cd )f(e cd ) J x f{e* d ) 

3 

=0(5)(1 + t) _ 2 e _ T+r, as x + 00 , 


jj (0 cd )dx J - pd x d t [In (f(d cd ))] 


and 


(f'(0 cd )d t 6 cd ) 2 cixf 

n 2 = -p- — f(Q C dj — = e 1+t > as x +°°> 


TZ 3 = d x ( ) = 0(<5)(1 + t) 1 e i+t , as i-a + 00 . 


1.3 Main results 

Now we are in a position to state our main results. 

Theorem 1.1. For any given [u+, u+, 0+,p_] with iq_ > 0 and 9+ > 0, suppose that [v-,U-,6-\ satisfies 
(11.181) . <j>± = pf 1 (v±) with <j>± £ and the function p e (-) satisfies the assumption (A). Let 

\v cd , u cd , 9 cd , </> crf ] (x, t) be the viscous contact wave defined in (|1.18D and (|1.21D with strength 5 = |0+—0_|. 
There exist positive constants eo > 0 and Cq > 0, such that if [i'o(x) — v cd (x,0),uo(x) — u cd (x, 0)] £ H 1 , 
[0o(aO — 9 cd (x, 0)] £ Hq and 

|| [u 0 (x) - v cd {x, 0), u 0 (x) - u cd (x , 0),9q(x) - 9 cd (x, 0)] || ffl + <5 < e 0 , 

then the initial boundary value problem <dZD> m and dUD admits a unique global solution [x, u, 9 , if] (x, t ) 
satisfying [x — v cd , u — u cd \ £ (7(0, + 00 ; H 1 ), \0{x) — 9 cd , (j> — (f> cd \ £ (7(0, + 00 ; Hq) and 

sup || [v - v cd , u - u cd , 9 - 9 cd , - qb cd } || < C 0 eo /3 - 

£>0 

Moreover, it holds that 

lim sup | \v - v cd , u - u cd , 9 - 9 cd , (f - <f cd 11 = 0. 

t-H-oo xgR+ 


( 1 . 22 ) 

(1.23) 
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From a physical point of view, the motion of the ion-dust plasma (cf. [321 US]); the self-gravitational 
viscous gaseous stars (cf. [3] ) and the charged particles in semiconductor devices (cf. [37]) can be governed 
by the NSP system. On the other hand, the NSP system at the fluid level can be justified by taking 
the hydrodynamical limit of the Vlasov-type Boltzmann equation by the Chapman-Enskog expansion, 
cf. n [nun ns]. In recent years, there have been a great number of mathematical studies of the NSP 
system. In what follows, we only mention some of them related to our interest. Ducomet E3 obtained 
the existence of nontrivial stationary solutions with compact support and proved the dynamical stability 
related to a free-boundary value problem for the three-dimensional NSP system in the case that the 
background profile is vacuum. Donatelli [5] established the global existence of weak solutions to the 
Cauchy problem with large initial data. Recently, Ding-Wen-Yao-Zhu [7] proved the global existence 
of weak solutions to the one dimensional isentropic NSP system with density-dependent viscosity and 
free boundary. Donatelli-Marcati [9] studied the quasineutral limit by using some dispersive estimates of 
Strichartz type. We point out that some nonexistence result of global weak solutions was also obtained 
in Chae [2]. Zhang-Fang m studied the large-time behavior of the spherically symmetric NSP system 
with degenerate viscosity coefficients and with vacuum in three dimensions. Jang-Tice [281 investigated 
the linear and nonlinear dynamical instability for the Lane-Emden solutions of the NSP system in three 
dimensions under some condition on the adiabatic exponent. Tan-Yang-Zhao-Zou [44] established the 
global strong solution to the one-dimensional non-isentropic NSP system with large data for density- 
dependent viscosity. In the case when the background profile is strictly positive, the global existence 
and convergence rates for the three-dimensional NSP system around a non-vacuum constant state were 
studied by Li-Matsumura-Zhang [33], Zhang-Li-Zhu [46] and Hsiao-Li [51] through carrying out the 
spectrum analysis. We point out that Duan m also used the method of Green’s function to obtain the 
large time behaviors of the more complex Navier-Stokes-Maxwell system. 

Another interesting and challenging problem is to study the stability of the NSP system on half 
space, to the best of our knowledge, there are very few results in this line. Duan-Yang m recently 
proved the stability of rarefaction wave and boundary layer for outflow problem on the two-fluid NSP 
system. The convergence rate of corresponding solutions toward the stationary solution was obtained in 
Zhou-Li [35| . We remark that due to the techniques of the proof, it was assumed in [13] that all physical 
parameters in the model must be unit, which is obviously impractical since ions and electrons generally 
have different masses and temperatures. One important point used in m is that the large-time behavior 
of the electric potential is trivial and hence the two fluids indeed have the same asymptotic profiles which 
are constructed from the Navier-Stokes equations without any force instead of the quasineutral system. 
Duan-Liu m then improved the results of [13] in the sense that all physical constants appearing in the 
model can be taken in a general way, and the large-time profile of the electric potential is nontrivial on 
the basis of the quasineutral assumption. For the investigations in the stability of the rarefaction wave of 
the related models, see also [12] for the study of the more complicated Vlasov-Poisson-Boltzmann system 
with more general background profile. 

When there is no self-consistent force, the NSP system reduces to the well-known Navier-Stokes 
equations. It is known that there have been extensive investigations on the stability of wave patterns, 
namely, shock wave, rarefaction wave, contact discontinuity and their compositions, in the context of 
gas dynamical equations and related kinetic equations. Among them, we only mention [H [23 [23 EH 
[Mi nil [an [Hi imi eh uni hh [mi [m l45l and reference therein. Moreover, we would also point out 
some previous works only related to the current work. Huang-Mastumura-Shi ]24| proved the stability of 
contact discontinuity of compressible Navier-Stokes equations with free boundary for the ideal polytropic 
gas through the construction of viscous contact wave profiles, the key observation in [24] is that the 
asymptotic profile of the temperature 9 satisfies a nonlinear diffusion equation, which can be solved by 
the technique developed in [I] [22] j and later on Huang-Mastumura-Xin [25] and Huang-Li-Mastumura 
[23] established the stability of the contact waves of the Cauchy problem. Recently Huang-Wang-Zhai 
[26] extended the results in [2!] to the general gas, however, for the Cauchy problem, it still remains an 
interesting open problem to generalize the results in [25j [23] for the general gas. 

In this paper, we intend to study the stability of the contact wave of the NSP system m with 
free boundary. Motivated by m and [21], we first construct the nontrivial asymptotic profiles of the 
quasineural Euler equations, it should be noted that the background density p e {4 >) satisfying assumption 
(4) allows that the asymptotic profile of the electrical potential can be distinct at the boundary. Then we 
perform the elementary energy estimates to the perturbative equations to obtain the global existence and 
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the large time behaviors. Compared to the classical Navier-Stokes system without any force, the main 
difficulty in the proof for the NSP system is to treat the estimates on the terms caused by the potential 
function (f). Precisely, the delicate term — d *td ) ip can not be directly controlled, as in El. the 
key point to overcome the difficulty is to use the good dissipative property from the Poisson equation by 
expanding p e (4>) around the asymptotic profile up to the third-order. In addition, it is shown [TT] that 
the sign of the first derivative of the rarefaction profile of the velocity and the good time decay properties 
of the smooth rarefaction profiles are important to the a priori estimate. Thus compared with m in 
which the stability of the rarefaction wave of the NSP system is proved, a new difficulty will arise, that 
is, the critical term (d x d cd ) 2 dxdt is beyond control, unlike that of [23], we need to pay extra 

effort to take care of the terms involving the self-consistent force, and it can be seen that the assumption 
(A 2 ) plays an essential role to obtain the desired estimates, see Lemma 1-01 for the details. 

The rest of the paper is arranged as follows. In the main part Section 2, we give the a priori estimates 
on the solutions of the perturbative equations. The proof of Theorem 11.11 is concluded in Section 3. In 
the Appendix, we present the details that are left in the proofs of the previous sections for completeness 
of the paper. 

Notations. Throughout this paper, we denote a generally large constant by C, which may vary from 
line to line. For two quantities a and 6 , a ~ b means < b < Ca. L p = LP(R + ) (1 ^ p ^ 00) 
denotes the usual Lebesgue space on R + = [0, + 00 ] with its norm || • II lp, and for convenient, we write 
|| • || L 2 = || • ||. We also use H k (k > 0) to denote the usual Sobolev space with respect to x variable on 
R + . C([0, T]; H k )[k > 0) denotes the space of the continuous functions on the interval [0,T] with values 
in H k . We use (-, •) to denote the inner product over the Hilbert space L 2 . [/1, / 2 ] £ H 1 means /1 £ H 1 
and / 2 £ H 1 , and so on so forth. 


2 The a priori estimates 


In order to study the stability of contact wave of the initial boundary value problem giD , m and 
(USD, that is, to prove Theorem ll.il we first define the perturbation as 

YP,^,Co-]{x,t) = [v - v cd ,u - u cd ,9 - 9 cd ,4> - cj) cd ] (; x,t ). 

Then [ip,i/j,(,a](x,t) satisfies 


d t ip - d x ip = 0 , 
d t i!> + d x p - d x p cd = 
R 


d x b ft: 


icd 


V V' 

cd o „, cd 


cd 


■ pd x 


d x i/) 

V 


7-1 

/ Ux 


~ x . r, car. ca a / <9 x <9 d x 0 cd 

dtC + pd x u - p d x u c = Kd x ( —- - 


V 


~F, 
+ G, 


dx ^ = -<P + V [1 - v cd p e (cr + (j ) cd )] - v cd d x i ^ x( ^ j , 

p(v , 9) - (0, t)=p-, C(0, t) = cr( 0 , t) = cr(+oo, t ) = 0 , 

b, C](^. 0) = [<£> 0 , to. Co] ( X) 

= [w 0 (ie) - v cd (x , 0 ), u 0 (x) - u cd (x , 0 ), 6>o(x) - 9 cd (x, 0 )] 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 

( 2 . 6 ) 


where x > 0, t > 0, F = —dtu cd + pd x ( dx “ ) and G = We note that the structural identity 

(12.411 will be of extremal importance for the later proof. 

The local existence of (ED. ED and ED can be established by the standard iteration argument 
cf. |21] and hence will be skipped in the paper. To obtain the global existence part of Theorem 11.11 it 
suffices to prove the following a priori estimates. For results in this direction, we have 

Proposition 2 . 1 . Assume all the conditions listed in Theorem 11.11 hold. Let [</>,?/>,£, tx] be a solution 
to the initial boundary value problem ED. ED- ED- (ED- (ED and ED on 0 < t < T for some 
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positive constant T. There are constants S > 0, eo > 0 and C > 0, such that if [ip,if] £ C(0,T; H 1 ), 
[C,cr] £ C(0,T; Hq) and 


sup ||[^,'0>C.o'](t)|| ff i +6 < e 0 , 

0 <t<T 

then the solution [ip,if,f,a](x,t) satisfies 

sup \\[T,ip,(,v}(t)\\Hi + [ \\dxV \\ 2 + \\dx [if,(, c] llffi dt 

0 <t<T Jo 

<CS + C ||[y>o> V’Oj Co]||^i 3 ■ 


(2.7) 


( 2 . 8 ) 


Proof. We divide it by the following three steps. 

Step 1. The zero-order energy estimates. 

Multiplying (|2.ip . (|2 . 2p and (12.31) by —R9 cd (L — ^), if and (9 1 , respectively, then taking the 
summation of the resulting equations, we obtain 



h 


where 


<h(s) = s — 1 — Ins, 

H = (p - p cd ) if - fj.- 

u 

Q ■ = (£) ~ ( 2 - fs ~ v) + (t) + ? ^ 


„cd\ afe r ^( dx ° 9x9 

v 9\v v cc 


and 


d x 9 _ . <pd x C Q , icd 






. C^x 6 » n d 


9 2 vv cd 


d x 9 ca . 


Let us now consider the most delicate term I\ on the right hand side of (12.91) . The key technique to 
handle I\ is to use the good dissipative property of the Poisson equation by expanding p e (<r + (f cd ) around 
the asymptotic profile up to the third-order. Only in this way, we can observe some new cancelations 
and obtain the higher order nonlinear terms. 

With the aid of (12.41) and (12.11) . one has 


_ d x ifa ifd x (pa tfd x v cd a ifd x v cd ip 

il —-1-o-1-o-r 


d* I — 


= -d. 


(v cd ) 3 v p' e (<f cd ) 

-V cd d x + v (1 - v cd p e (a + cf cd )) - v cd d x ( 


d x <f‘ 


cd\~ 


h,i 


J-d x 


-v cdd x + V (1 - v cd p e (a + cf cd )) - v cd d x ( 


d x (f‘ 


cd\ 


ifav 


1 1,2 

- v cd dX ( + V (1 - v cd p e (a + cf cd )) - v cd d x ) fidxV^r^r)]- 1 


( 2 . 10 ) 


ifd x v cd a f a'lf 
- o -+9 ; 


12 


V 





























To deal with the lower order terms involving 1 — v cd p e (cr + (j ) cd ), we first get from the Taylor’s formula 
with an integral remainder that 


1 - v cd p e {a + <t > cd ) = - v cd p' e {<j> cd )o - 


v cd p'M d ) 0 „. cd 


Pe'(Q) 


(</> - Q? 


dg. 


Io 


By virtue of (12.111) . we then compute /ip, /12 and /13 as follows: 


J M=-^ 


+ 


,cd 


-(9 x cr ) 2 - dt 


d x 


v cd \ d x a 


d t a - v cd d, 


-ad x 

d x cr 

V 


,cd 


~2 dt 


.cd 


( d x a ) 2 


ad t (v- 1 )+-d t (v cd p'M d )a 2 ) 


+ \dt (v^p'^a 3 ) - dt I 0 a+^v cd p'M d )a 2 + \d t (v cd p' e {y cd )) a 2 
6 v 2 


I 3 


1 d t v yCd Jt( A,cd \^3 

2 V 


+ d x d t 


p"(<TV + ^ (t^p"^)) a 3 - I 0 a + d t (v cd d : ( ^ 

)-«■( 


(TV 


-1 


fv cd ad x cr 

V 7 2 


( v cd d t crd x a 


..cd 


..cd 


h ,2 =-^-d 2 ad x ay + — d 2 aad x (tpv 2 ) + v cd d x crd x (v L )d x (a'tj;v 2 ) 

..cd j , d j ,icd ,d x vcr 2 il> r.cdjficd^v 2 ^ „.cdjuicd^ d x ( r <j 2 i } 


-v cd p'M d )- 


-v cd p'M d )- 


- d x (v cd P 'M d )) — -v cd P 'M d )- 


> 2 V V v z 


2 

+ v cd d, 


d x fi 


cd 


h,3 =~d, 


+ 


d x ( y<rv 2 ) - d x (v cd d x y av 2 ^j > 

yd x v cd d x v cd ay 1 p'l{4> cd )d x v cd a 2 y 


d x 

v J v(v cd ) 2 p' e ((f> cd ) ( v cd ) 2 2 (y cd ) 2 p' e {(j) cd ) 

d x v r I 0 ip 


h 


I e 


(v cd ) 3 p' e ((j) cd ) 


- d x 


d x (j)' 


cd 


ipd x v‘ 


cd 


( 2 . 11 ) 


( 2 . 12 ) 


(2.13) 


(2.14) 


v J v(v cd ) 2 p' e (cj) cd ) 

Note that /; (2 < l < 6 ) can not be directly controlled. To overcome this difficulty, we first get from (12.41) 
and ( 12 . 111 ) that 


{I2 + I5 ) + 1% + Z4 — — 


d x v cd a4np(v + v cd ) p"{<j) cd )d x v cd <j 2 y cd , d d x a cry 


( vv cd ) 2 


2 p' e (4> cd )(v cd ) 2 


„ca f / ica\ 

-V p e {y )- 




V* 

d x v cd a 2 y 


(y cd ) 2 

d x va 2 y n CdJ (xcd ^ x 1 d x v cd a 3 y(v + v cd ) JlfACd ^ d x v cd ai/jl 0 {v + v cd ) 
v 2 

d x v cd a 




d x v cd a 2 y p''{y cd ) 

V V cd p' e {(j) cd ) 


v ;(0 + o 


vv 


cd 


p'M d ) - 


+ 


yd x ( 


v(v cd ) 2 


dx<t> c 


v + v c ) ' d xV cd ayd x (^f-) (v + v cd ) 

I o • 
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which is further equal to 


1 d x i(>(T 2 

2 v 


v cd p'M d )+^v cd p'M d )d x - -7 +o 


1 1 


V V 


1 a 2 tl>ip 


2 v 2 


d x v 


cd 


p'M cd ) - 


VPM^P'L^) 

V cd p' e {(j) cd ) 


17 


l a^dxp ^ d 1 d x v cd a 3 ip(v + v cd ) cd _ d x v cd ^I 0 {v + v cd ) 
2 v 2 1 PeW ) ' 2 ~.~.rA ^eVS° / 


VV { 


v(y cd ) 2 


(2.15) 


d x v cd atpd x {v + v cd ) d x v cd a^d x (& 2 ) (v + u cd ) fi/> 


For J 3 and I 7 , it follows from (Or, d and (12.41) that 


h + h = 


2 d T u cd a 2 


Sd Wr „ dJ ,, cd 


2 v 


Veit'*) 


=d t 


3pa 2 

2 v 


~ d *[ I ^v cd p' e (cj> cd ) 


■ cd P 'M> cd ) 


V Cd p'M Cd ) - S d p'e 


(2.16) 


o 2 f) n ,cd—2 

^ o ( n ,°d rJ (A>. c d\\ u cd / / icd\ 

~ 2 dt \ V Pe\9 )j — + --- v PeW )• 


Plugging (12.161) . (12.151) . (12.141) . (12.161) . (12.121) and (12.101) into (12.91) . integrating the resulting identity with 
respect to x over R+, and using (A 2 ), we thus arrive at 


d_ 

dt 




(M, 

TD 

pcd*i 

VtW ^ 

7-1 

w cdN 

3 d f 

- ) da; - 

- / 

* y 

^ dt J r + 


e cd 


cd cd 

+ ~ 2 ~\p'M cd )\a 2 + ^{d x a) 2 )dx 


+ p f dx ^ dx+ f d x () 2 dx 

Jwl+ v J R+ v 6 

31 

— f Qidx — ( Q 2 dx + f Fipdx+ f ^Gdx + H( 0, t) + 

J R_l J Rx J R , «/R_i_ P 1 1 


(2.17) 


where 


H = H(x,t) ={p - p cd )il> - p 

v cd d t ad x a 


d x ^ C 

-ft~7 


3*0 

V 


d x ® 


cd 


nCd 


- 8 , 

V 


( v cd ad x <x 

V ^ 




Ml + 2 + 


V ) 


2 v 


v cd P 'M d ), 
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and 


h = / d x 


v cd \ d x a /' cd / d x cr\ 

dtcrdx , X 2 = — / v a x ( - ) ox/t (v )ax, 


v J v 


+ 
cd > 


2 3 = - / dJocrdx, I 4 = - 


z 5 = - 
2 7 = 

29 = 
2 n = 
Zl3 = 


d t v 


d x r d 


I 0 adx, 1 6 = dt[ v ca d, 


'j 1 dx, 


f cy f 9 

/ —o-d x ad x <T'tpdx, Xs= — d 2 aad x {ypv~ 2 )dx, 

JR+ v JR + v 

[ v cd d x ad x (v~ 1 )d x (<j'i/jv~ 2 )dx, Iw = \ [ d t {v cd p' e ((j> cd ))a 2 dx, 
JR+ Z J R+ 

1 f ^L v - d p'^ d )aHx, 1 V 2 = [ 

2 Jr, v Jr. 


_d x y^ d I^_ 

(v cd ) 3 p' e {(f> cd ) ’ 


d x vl 0 cnl 1 


dx, I 14 = 


dxlo&ip 


dx, 


X 15 — — 


Ii7 = - d : 


d x 


il>d x v' 


cd 


V J v(v cd ) 2 p' e (4> cd ) 
d x (jf d \ ipd x v cd 


dx, lie = / u cd <9 x 


d x r d 


d x (%/jcrv 2 ) dx, 


+ \ V J v(v cd ) 2 p' e (4> cd ) 

1 /' ..cdn^cd^VCF 3 ^ 


dx, 2 18 = - / 5 t (v cd p"{(j) cd )) a 3 dx, 

o ./TO , 


T _ 1 /' ..cdJ'fJ.cd&Vril) t 1 /' a (..cd J't J.cd\\ V 

-^19 o / ^ Pe \4* ) o -^20 o / ^ P e 

z JR + u ^ JR + v 

[ v - d p" e {^ d ) d J^± dx , x 22 = [ a 2 ^v cd p'M d )d* (~ - -4) <**> 

JR + w JR+ \ w w / 


221 =- 


1 /' „,cd 


2 23 ~ ~ 


i«=-u 

2 7r + ^ 


/ / ,cdN _ ^Pe(0 Cd )Pe(0 ) 

6 V cd p' e ((j) cd ) 

p' e {(p Cd )dx, l25 = \ [ 
z «/R 


dx, 


1 /• + O d 

— PeVP ) d X, 


R+ 


,,cd 


^26 = — 

2 28 = / 

Ji 

230 = - 


d x v cd ai>Io{v + v cd ) 


R + v(v cd ) 2 

d x v cd ai/jd x (&f-) (u + u cd ) 


dx, 2 27 = 


d x v cd a^d x (u + u cd ) 


dx. 


dx, 2 2g = 


d T u cd " 2 


v cd p'M cd )dx, 


1 / a, (v cd p'M d )) — dx, 2 31 = -|/ v cd p' e {<j) cd )pd t (—') dx. 
2 Jr, ^ 2 J r V v / 


We now turn to estimate the right hand side of (12.17|) term by term. It should be noted that the following 
Poincare type inequalities play an important role in our computations: 


I COM) I < £ 2 ||0*CII, Iv’Ca:, i)| < |<p(0,t)| +x“ ||5 X ^||, |cr(x,t)| < x= ||a x cr|| 

From (12.181) and Lemma 14.11 one can further obtain 

f P 2 ((d x 0 cd ) 2 + |a 2 d cd |) dx < CS 2 Wpofnie-^* + CS 2 \\d x pf, 

J R_l_ 

/ (C 2 +a 2 ) ({d x e cd f + \d 2 x e cd \) dx < CS 2 \\d x [C,a}\\ 2 , 

JR + 

where the following Sobolev inequality is also used: 

|d(x)| < •\/2||/i|| 1 ^ 2 ||9 x d || 1 / 2 for h(x) £ H l ( K + ). 


(2.18) 


(2.19) 


( 2 . 20 ) 
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By applying (12.191) . Lemma f4.21 the a priori assumption (12.71) . Cauchy-Schwarz’s inequality with 0 < 77 < 1 
and Sobolev’s inequality ( 12 . 201 ) . we obtain the estimates for terms involving Q\ and Q 2 as follows: 


Qidx 


<C f (if 2 + C 2 ) ((d x 6 cd ) 2 + \d 2 0 cd \) dx < + G(5||c0[<^,C]|| 2 , ( 2 . 21 ) 

Jr 


Q 2 dx 


<(Ce 0 + V )\\d x C \\ 2 +0^1 (<p 2 + ( 2 )(d x 0 cd ) 2 dx 
J R + 

<(Ce 0 + r))\\d x (\\ 2 + C^WipoW^e^ — * + C^WdxlipXW ■ 


( 2 . 22 ) 


For the terms involving F and G, noticing that 

\d t u cd \ = 0(l)tf(l + \d 2 x u cd \ = \d x d t v cd \ = 0(1)<5(1 + , as x -> + 00 , 

we get from Cauchy-Schwarz’s inequality that 


< f \d t u cd ip\dx + C f \d 2 u cd ijj\dx + C f \d x u cd d x v cd tp\ dx + C f \d x u cd d x <pi/j\ dx 
Jr+ J k _(_ Jr+ Jr + 

<CS(1 + t)~ l ~ a ||0|| 2 + C8{ 1 + 0~ 3/2+ “ + C5\\d x ip\\ 2 , 

(2.23) 

where 0 < a < 1 / 2 , and 



'-Gdx 


< G||C|U|cU || 2 < C'eoll^V’ll 2 + Ce 0 (1 + 1 )~ 


(2.24) 


We next compute the term H(0,t) arising from the boundary. Since C(0,t) = cr(0, f) 
be reduced to 


0 , 




can 


R6-<p{0,t) 

v( 0 ,t)v- 


■ 0(0 ,t) + /u, 



(o,0 


which is further dominated by 


G|v?(0, t)0(0,01 + C|(<9 t v0(0, 00 ( 0 ,01 

<^1^0(0)16-^*10(0,01 < G||^o(^IUOl0ir / 2 |l^0l| 1 / 2 e"^ t (2.25) 

<G||^o(x)||^ / 1 3 e-^‘+Ge 2 ||a x 0|| 2 , 

according to Lemma 14.21 Sobolev’s inequality (12.201) and Young’s inequaity. 

In order to estimate Ii (1 < l < 31), we first calculate 


Jo ~ a 3 , dtlo = -v cd d t <p [* (0 - e) P :\e)d B + \cr 2 d t <f d v cd P ":^ d ) - dtv cd r (0 o eY P "(o)dQ 

J (fi c d " J (j) c d " 

~ 9 4 0ct 2 + d t v cd o 2 + d t v cd <j 3 = d t oa 2 + 2 d t v cd cr 2 + d t v cd a 3 , (2.26) 

and similarly, 

d x I 0 ~ d x cra 2 + 2 d x v cd a 2 + d x v cd a 3 . (2.27) 

In addition, from <j2.4|) and (12.5[) . it follows 

\\d t a \\ 2 + \\d t d x a \\ 2 < C||0 B 0|| 2 + Ce 0 |j [d x <p, 5 2 0, d x a, d 2 a] f + CS(1 + t)~i . (2.28) 

For the sake of completeness, the proof of (12.281) is given in the appendix. 
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With (12.261) . (12.271) and (12.281) in hand, we now employ (12.191) . Cauchy-Schwarz’s inequality with 
0 <r] < 1, Sobolev’s inequality and Lemma ITT! repeatedly to present the following estimates: 


12d | \d x [v cd ,v\d x crd t a\ dx < Ce 0 || [ d x a,dla,d t cr] ||“ , 

Jr + 

\X 2 \ < C f \dlaad x [4>,u cd } \ dx + C f \d x aad x [ip,v cd ]d x [ip,u cd }\ dx < Ce 0 \\d x [(p,ip,a,d x 
Jr+ Jr-l. 


|X 3 | + |2 5 | + |Xn| + |Iis| 


<C f \d x ii cd a 3 \ dx + C f \a 3 [d t a, d x ip]\dx + C f \d x ipcr A \dx + C f \d x 
JR+ JR+ JR-\. Jr 

<Ce 0 \\[d t (J,d x a, d x ip ]\\ 2 , 

|J 4 | < C [ \d t {v cd ,v]\\d x a\ 2 dx < Ce 0 116*^[cr, cr]11 2 , 

JR+ 

'd x r d ' 


i cd cr A I dx 


|2e|< 


d t d x v 


av 1 

dx + 

Q Cd ( d x4 )Cd \ o -1 

d t v - d x av 


J R+ 

\ v J 


dx 


led 


d t v c " ( —— ) ad x (v 1 ! 


+ 


..cdfl I'dxfi 


v cd dt 


cd 


d x av 1 


dx + 
dx + 


d x v cd d t 

„ ft 
v cd d t 


d x <t> { 


cd 


cd a 1 9 x r d \ , 


ad x (v ) 


dx 

dx 


<C 6 (l+t) ^ + Ceo\\d x [ip,4>,d x <r}\\ 2 , 

12 | ||p2„||2\ 


|I 7 | < c\m L ~(\\d x *r + IIW) < Cca\\d x [(j, d x a]\\ 2 , 

1^81 <C* f \d 2 crad x tp\dx + C f \d x aaipd x ip\dx + C f \d 2 aaipd x v cd \dx 

J K_|_ J M_(_ J M_|_ 

<Ce 0 \\d x [ip,ip,d x a}\\ 2 + Ce 0 f (d x 9 cd ) 2 cr 2 dx < Ce 0 \\d x [ip, tp, a, d x cr]|| 2 , 

JR+ 

\Xq\ <C / \d x crd x vdxcr'il)\ dx H- C / \d x (jd x vcrd x f il)\ dx + C / | d x ad x va'ifjd x v | dx 

J R + J R + Jr + 

<C\\d x a\\m IMIffi lino’ll \\d x v\\ + C\\d x a\\ H i \\cr\\ H i \\d x ip\\\\d x v\\ + C\\d x a\\ H i ||cr|| ff i \\ip\\m ||3x<p|| 2 

\a\\ H i\\^\\ H i\\d x v cd \\ H i+C\\ip\\ H i (||9 x a|| 2 + ||(7(9 x i; cd ) 2 || 2 ) 


+ C'||ct c ij|| 

<Ce 0 \\d x [a,d x (j, ip,ip}\\ 2 , 


|Xio| + | 229 1 + |2 30 | < C f \d x u cd a 2 \dx < C f (( d x 9 cd f + \d 2 J cd \) a 2 dx < C5\\d x a\\ 2 , 
Jr+ Jr+ 


|2i 2 | + |2i 3 | + |Zi4| + |2ig| + |X2o| + 1^25 | + 1^261 


<C f \d x v cd a 3 ip\ dx + C f \d x [ip,a]a 3 tp\ dx 
J r + J R + 

<Ce 0 \\d x [a, ¥>]|| 2 +Ce 0 [ (d x 6 cd ) 2 a 2 dx < Ce 0 \\d x [a,tp]\\ 2 , 

J R + 

| 2 i 5 | <rj f \d 2 a\~ dx + C v ( ip 2 (d x 6 cd ) 2 dx + C f \d x crd x tp\ 2 dx + C f \d x ad x v cd \ 2 dx 

Jr -|_ J M_)_ Jr.\- Jr -|_ 

<{Ce 0 +r))\\d x [a,d x a\\\ 2 + Cr, f i) 2 {d x e cd ) 2 dx, 

J M_l 
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|Xi61 <C f \d 2 v cd d x aip\ dx + C f \d^.v cd ad x , tp\ dx 

J M_l_ J R_|_ 

+ C f \dlv cd [d x tp,d x v cd }cnp\ dx + C f \d x v cd [d x f,d x v cd ]d x aip\ dx 
Jr+ JR + 

+ C f \d x v cd [d x p,d x v cd ]ad x ip\ dx + C j \d x v cd [d x ip, d x v cd ] 2 aip\ dx 
J M_|_ J M_|_ 

<Ce 0 \\d x [(j, if,ip ]|| 2 +C5{l + t)~i, 

|2 17 | + |2 27 | <C / \{d x v cd ) 2 dM\ dx + C [ \d 2 v cd d x v cd ^\ dx 
j\ r_l Jwu 


<Ce 0 \\d x p\\ 2 + C / i/j 2 (d x 0 cd ) 2 dx + C 6 2 (l + t)— 2 , 

J K + 

12211 + |2 22 | + |X 23 | + 12241 <C [ \d x [a,if]a 2 ip\dx + C f \a 2 ipipd x v cd \ dx 

JM.+ */R_|_ 

^Ceoll^Io-, ^]|| 2 +Ce 0 f {d x 6 cd ) 2 ip 2 dx 
Jr+ 

<Ce 0 ||a x [cr, ^]|| 2 + C5 2 \\f 0 \\ 2 H ie~~ t , 

|2 2 81 <C [ \d 2 ad x v cd aip\ dx + C f \d x ad x vd x v cd ipa\ dx 
Jr+ Jr + 

<Ce 0 \\d x [a,d x a, p }\\ 2 + Ce 0 f (d x 0 cd ) 2 a 2 dx < Ce 0 \\d x [a, d x a, p]\\ 2 . 

J M_l_ 

For the last term Z 31 , applying (12.111) again, one can see that 


2 31 =- 


\ jf v cd (JM> cd )d t ^ j (-v cd d x + t; (1 - v cd p e (a + <j > cd )) - v cd d x ) 


dx, 


which implies 


2 31 - 


dt 


[ (v cd p' e {<t > cd )) 2 a 3 dx 
J M-L 


< Ce 0 || [ d x ^,d x f,d t (J,d x a,d 2 a ] ||~ + C5{1 + t) = . 


Let us now define ao(x) = cr(x, 0) = (p(x,0) — (p cd (x,0). From the Poisson equation (12. 41) . it follows that 
for any t > 0 


\Wmh<cMt )\\ 2 + c\\d 2 x ^ d ( 

and hence in particular, 


+ C 


(■ ®A cd ) (t) <C\\p(t )\\ 2 + C5 2 , 


(2.29) 


\Wo\\h <c\\p 0 \\ 2 + c5 2 . 

We now conclude from (12.171) . (12.211) . (12.221) . (12.231) . (|2.24|) . (1 2.2 5 [) . (12.281) . (12.291) and the above 
estimates on I; (1<Z<31) that 


m,<P ,{]\\ 2 + M 2 m +e 0 \\d x <p \\ 2 


1 ll&MII 2 * 

0 

<C||[V , o,Co ]|| 2 + C'||^o||^i 3 + (Ce 0 + 7 ?) [ \\d x [p,d x -ip,(j,d x <j}\\ 2 dt 

Jo 

+ C V [ [ ip 2 (d x Q cd ) 2 dxdt + CS, 

Jo J tx 


(2.30) 


for suitably small eo > 0, S > 0 and p > 0. 
Step 2. Dissipation of d x [ip, a, d x a]. 
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We first differentiate m with respect to x, to obtain 


d t d x p - dl^) = 0. 


(2.31) 


Then multiplying (12.41) . (12.21) and (12.311) by c^er, —vd x <p and pd x tp, respectively, and integrating the 
resulting equalities with respect to x over R + , one has 

/ -— {dl(i) 2 dx+ f vv cd \p' e ((j) cd )\(d x a) 2 dx 

J R+ v J R+ 

= [ ^-d x vd x adladx - f d x v [l - v cd p e (cr + (f> cd )\ d x adx + [ vd x [v cd p' e (<j> cd )]ad x (jdx 

Jr+ v J R+ J R+ ( 2 . 32 ) 

d x adx + / d x pd x adx + ip(0,t)d x a(0,t) — / vd x Iod x adx 

J Rj. J Rx 


+ / vd x 

J E_l_ 

..cd 


v cd P ':(r d ) 2 


V cd d x 


d x (t> [ 


cd 


dladx , 


cd j.cd 


d x <j) cd d. 


v v l 


/ d t il>vd x tpdx - / (d x p - d x p cd )vd x ipdx + / d x tpd x crdx + / 

J R+ Jr + J R + J R_! 

= — / pdlipd x pdx — / pdlu cd d x tpdx — p / d x (v~ 1 )d x uvd x ipdx + / dtu cd vd x ipdx , 

«/]R_i_ 


vd x pdx 


and 


/ /r (d t d x tp - dli/)) d x pdx = 0 . 
J R_l_ 

The summation of (12.321) . (12.331) and (12.341) further implies 

tl>vd x pdx + 7 :-r [ (d x p) 2 dx+ f p cd (d x tp) 2 dx 
i J r + 2k + 

+ [ — — (dla) 2 dx + j vv cd \p' e {(j) cd )\ld x (j) 2 dx 

ilj. v J Ex 


(2.33) 


(2.34) 


d_ 

dt 


f Rd x 


vd x pdx — / Rtpd x 

■ fled - 

/R+ 

M 

J M l 

VV ca 


vd x tpdx 


J\ J 2 

d x (j) cd d x (j) cd 


J 3 


J 4 


V V 


cd 


vd x <pdx — / fid 2 u cd d x <pdx — p / d x (v 1 )d x uvd x ipdx 
JR+ ila 


J 5 






(2.35) 


» mCd 


+ 


d x vd x ad 2 adx — f d x v [l — v cd p e (a + </> cd )] d x adx + 99 ( 0 , £)<9;cCr(0, £) 

./r. v ---' 


J 8 


J 9 


[ vd x I 0 d x adx+ f vd x \v cd p' e {(j) cd )\ crd x adx + [ vd x 


j 10 

v cd P "M d ) 2 


d x adx 


J 11 


J12 


^13 


+ / d t u vd x pdx- / s i 

dR+ JR + 


c^ cd 


dzadx. 


J 14 


^15 


We now turn to compute </; (1 < l < 15) term by term. For brevity, we directly give the following 
computations: 

|Ti| <C f \i/)d x ipd x ip\ dx + C I \^d x u cd d x ip\dx<Ceo\\d x \i/),ip]\\ 2 + C6{l + t)~ 3/2 , 

J r_l Jrx 
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IJ 2 I < \ip(0,t)v(0,t)(dtip)(0,t)\ + 

[ v(d x ip) 2 dx 

+ 

/ ipd x vd x ipdx 




J ]R_|_ 


<C\\vo\\m e " + Ge o\\d x [ip,<p\\\ + C\\d x ip\\ +C ip (d x O c ) dx , 

| J 3 I + \ Ja\ + |«^ 5 1 < {v + Ceo)||<9 x (p || 2 + C' ?) ||9 X CI | 2 + C 17 [ (jp 2 + C 2 ){d x 0 cd ) 2 dx 

J R + 

< (ry + Ceo)||9 x [v ? ;C]l| 2 + Cr^ll^CII 2 + ^5\\tpQ\\ 2 H ie » *, 

|J6|<CT||5 x ^|| 2 + (7<5(l + t)- 5 / 2 , 

|-^r| / (d x ip) 2 \d x ip\dx + C j (d x ip) 2 \d x u cd \dx + C [ \d x v cd d x u cd d x ip\dx + j \d x v cd d x ipd x ip\dx 

J M_(_ J ]R_|_ M_|_ J M_(_ 

<Ce 0 \\d x [ip,d x ip,ip]\\ 2 + C5(l + t)~ 2 , 

|Js| + |Jg| < Ce 0 \\d x [y,d x (j,o\\\ 2 , 

IJ 10 I < Ce 0 \\d x a\\ 2 Hl + C'||^o|||ie- i ^ t , 

| Jn| + | J 12 I + | Ji 3 | + | Ji 4 | + | Jis| < Ceolia*^, <r, d x a}\\ 2 + C5{ 1 + 1)~ 3 / 2 . 

Substituting the above estimations for J; (1 < l < 15) into (12.3511 . letting 77 > 0 be suitably small and 
combing (|2.30D . we obtain 


^^C ]|| 2 + IHIhi + IM 2 + / \\d x a\\ 2 m dt+ \\d x [ l p,ip,(]\\ 2 dt 

Jo Jo 

<Ce 0 / T ||^|| 2 dt + Cd + C||[i/;o, Co]|| 2 + C ||^ollt/f + C P [ ip 2 {d x O cd ) 2 dxdt. 
Jo Jo J Rj. 


(2.36) 


Step 3. Higher order energy estimates. 

Multiplying (12.21 ) by —d 2 ip, and integrating the resultant equality with respect to x over R + , one has 


~[ (d x ip) 2 dx + [i [ dx 

2 dt J^ + 7 r + v 


4* V 


Led 


V V 


cd 


d 2 ipdx + [ d x {p — p cd )d 2 ipdx +p f 
J ij. Js.. 


d x ipd x ip „ 2 1 1 

- 2 - d xWx 


J 16 

/* d ibd f 

+/i / ——J- d 2 ipdx- / Fdlipdx—(d t ipd t ip)(0,t). 

J R + u v --—v-' 


Jl7 


J 18 


(2.37) 


J21 


J 19 


J 20 


To obtain the estimates for J; (16 < l < 21), we use Cauchy-Schwarz’s inequality with 0 < 77 < 1, 
Sobolev’s inequlity ( 12 . 201 ) and (12.191) repeatedly to perform the calculations as follows: 


|Jie| <C f \d x v cd <pd 2 ip\ dx + C f \d x crd 2 ip\dx 
J K_l_ J M_|_ 

<(CS + rf)\\d x [ip,d x ip]\\ 2 + C v \\d x a\\ 2 + C'(5||</J 0 ||| f ie"~ i , 

\Jn\<C f \d x [(,<p\dlip\ dx + C [ | [C, y\d x {y, v cd ]d 2 ip\dx 
Jr + Jr + 

<(Ce 0 +ri)\\d x [ip,d x -ip}\\ 2 + C v \\d x [(, ^]|| 2 + CdH^ollffie - ”*, 
\M + \M<Ceo\\d x [iP,d x iP}\\ 2 , 
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|J 20 | <C [ \d 2 x u cd d 2 x ip\dx + C ( \d x u cd d x vd 2 i/j\dx + C [ \d t u cd d 2 ip\dx 

«/R_i_ J R_l 


<C5\\d x [(fi,d x ^]\\ 2 + C5(l + t) 2 . 
For the last term J 2 1 , in light of Lemma 14.21 we have 


d 2 i = -(dti>d t ip)(0,t) = -d t [(V>dtyO(M)] + y?o( 0 ) ^ " *, 

A* 


furthermore, it follows that 


|O9t<A>)(0,T)| < C^ 0 (0)V’(0,T)e " T < Ceoll^ollffie * J , 


and 


|(M^)(o,o)| < qv-o(0)^o(0)| < c {\\ m\ 2 h i + llwll 


ff 1 - 


(2.38) 


(2.39) 


(2.40) 


By virtue of (12.38|) . (12.391) and (12.40|) and carrying out the similar calculations as (12.25[l . we thereby 
obtain 


J2idt 


<Ce o / ||5^|| 2 dt + C'||^o||^ 1 +C'||^o||^ 1 +Ceo||^olk 1 . 


Plug the above estimations for J ; (16 < l < 21) into (12.371) . and recall (12.361) and (12.301) . then choose 
e 0 > 0 , (5 > 0 and 77 > 0 suitably small, to derive 


\\W,V,(}\\ 2 + \W\\m + \\dM\ 2 + \\dM' 


\\d x a\\ H idt+ \\d x [ip,ip,Q\\ dt+ \\d x ip\\ dt 


<C5 + C WCoW+ C\\Wo,<p 0 }\\ 4 I i? + C I I V> 2 (d x 0 cd ) 2 dxdt. 

Jo J R+ 

Similarly, multiplying (12.31) by —d 2 (,, and integrating the resulting equality over R+, we obtain 


(2.41) 


R d 

2(7 — 1 ) dt 


(d x () 2 dx + k f 

J Rj 


{d 2 x Cf 


dx 


= [ ( pd x u - p cd d x u cd )dlQdx + n f 
Jr + Jr , 

V d x6 cd \ d 2(- dx _ f Gd ^ dx 

J R-L 


^^d 2 ( dx + nf dxCd f d d 2 x (dx 


(2.42) 


n d x 
J R + 


where we have used boundary condition £(0,i) = 0. The right hand side of (12.421) can be handled as J; 
(16 < l < 21), the details of which we omit, therefore one can get from (12.421) and (12.411) that 


W,<P,C]\\h* + Mh*+ IM 2 m dt+ \\d x [<P,M\\ 2 dt+ \\d 2 x [^C)\\ 2 dt 

J 0 Jo Jo 

<C6 + C\\['po,il>o,<Zo]t 1 g + C [ f 2 {d x e cd ) 2 dxdt. 

Jo J R + 


(2.43) 


Finally, letting 5 > 0 small enough, combing (12.43[) and (14.31) in Lemma POl we obtain (12.81) as desired, 
this completes the proof of Proposition 12.11 

□ 
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3 Global existence and large time behavior 

We are now in a position to complete the proof of Theorem 1.1. 

Proof of Theorem \l.l\ In view of the energy estimates obtained in Proposition l2.il one sees that 

sup \\\p,ip,C,,a]{t)\\ 2 m < CS + C\\[4>oXo,Po]\\m ■ (3.1) 

n<r+<r "T v J 


Notice that 5 > 0 is a parameter independent of eo- By letting 6 > 0 be small enough, the global 
existence of the solution of the Cauchy problem (12.11) , (12.21) , (12.31) , (12.41) , (12.51) and (12.61) then follows from 
the standard continuation argument based on the local existence (cf. [21]) and the a priori estimate (12.81) . 
Moreover, m implies (11.221) . Our intention next is to prove the large time behavior as (11.231) . For this, 
we first justify the following limits: 


lirn \\d x [tp,-if,(\{t)\\ 2 L 2 = 0, 

£—> oo 


and 


lim ||<3 r (r(t)|| 2 = 0. 

£—>•+00 


To prove (I3.2[) and (13.31) . we get from (12.11) . (12.21) . (12.31) . (12.81) and (11.201) that 



d 

dt 


II + + + 


r+oo 


dt =2 


\{d t d x [ip,ip,Q,d x [p,ip,Q)\dt 


r + oo 

<C + C \\d x [p, if, C, 9 X [ip, c, o-]] || 2 dt < Too. 
Jo 


(3.2) 

(3.3) 


(3.4) 


On that other hand, (12.28[) . (12.361) and (12.81) yield 


/•+oo 

d O 

L 




\(d t d x <j,d x a)\dt < Too. 


(3.5) 


Consequently, (13.41) . (13.51) together with (12.81) gives (13.21) and (13.31) . Then (11.231) follows from (13.21) . (13.31) 
and Sobolev’s inequality (12.201) . This ends the proof of Theorem ll.il 

□ 


4 Appendix 

In this appendix, we will give some basic results used in the paper. The first lemma is borrowed from 

Ell- 

Lemma 4.1. Let d cd satisfy (11.181) . for \6 + — 9-\ = 5, it holds that 

f {d x e cd fdx < C6 4 (l + t)-%, [ {d 2 x e cd fdx < C5 2 {l + t)-i, 

J R + J R + 

[ (dl9 cd ) 2 dx <C5 2 (l + t)~i, [ x((d x 6 cd ) 2 + \d 2 x 0 cd \)dx<CS. 

J ix J ti. 


Next is the key observation from the boundary condition (12.51) . 
Lemma 4.2. It holds that 

T(0 ,t) = To(0)e _ ~ i . 

Proof. Since d x u cd (0,t) = 0, from (12.51) it follows that 

RQ- dt<p(0,t) 

■ P ~— , _, n ^ =P~, t> 0, 


(4.1) 


v- + ip(0,t) v- + p(0,t) 
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which implies 


dtp(0, t ) = -—ip(0, t). (4.2) 

(ED follows from (14.21) and the compatibility condition <^(0,0) = <^>o(0). This ends the proof of Lemma 

□ 

We now give the following estimates concerning the delicate term / / ( d x 0 cd ) 2 ip 2 dxdt. 

Jo J R + 

Lemma 4.3. Assume all the conditions listed in Provosition \2.1\ hold, then for any 0 < T < +oo, there 
exists an energy functional £{<p, if, C) with 

\£{<P,M\<C5 2 \\[tp,M\\ 2 , 


such that the following energy estimate holds 

£(<P,M(T)+ f [ (d x 6 cd ) 2 ip 2 dxdt<CS + C6\\g> 0 f H i+C6 f \\d x [<p,Tp,(,a,d x a\\ 2 dt. (4.3) 

Jo J R+ Jo 

Proof. Define 

w= f ( d y 0 cd ) 2 dy. 

Jo 

It is easy to check that 

||wM)||oo < CS 2 (l + t)~^, ||d t tu(-,i)||oo < C5 2 (l + f)“3. 

From (12.41) and (12.111) . it follows that 


<P 


1 


VV cd p' e {4> cd ) vv cd p' e (<p cd ) 


v cd d x 


8 x a\ , v cd p e {<f cd ) 


+ 


o 2 — Iq I v + v cd d ; 


d x <j>' 


cd 


M 


On the other hand, (I2.2|) can be rewritten as 

' R( - p cd p - a ' 


d t if + d x 


Substituting (14.51) into (14.61) . one has 


= ~d x - \<J - 


1\ tpd x (j>' 


cd 


VV 


cd 


pd x (^)+F. 


d t if + d x 


'R(+ -P cd )g> 


M 


= d x ( — ]-d x \-)a- 


1\ pd x (p' 


cd 


VV 


cd 


(4.4) 


(4.5) 


(4.6) 


+ pd x (^)+F. (4.7) 
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Multiplying (14.71) by RC + — P cd ) P vw , integrating the resulting equation over R + leads to 


RC + 


1 


d 

dt 






R( + 


RC ■+■ 


VV cd p' e ((f> cd ) 

1 


p cd l<p 


cd 


i p , M> cd ) 


i 


(. d x e cd ydx 


vwdx — / ipdt 

J R + 


RC 4" 


" cd p' e {(j) cd ) 


-p cd \p 


vwdx 


w cd p' e {(j) cd ) 


P cd \<P 


dtvwdx — / ip 

J Ro- 


RC + 


K 1 

1 


W cd p' e ((f> cd ) 


P Cd )<P 


vdtwdx 


K.2 


RC + 


n cdp / 


- P cd lP 


wdx +p [ -^-d x 

J R+ V 


K.3 

RC + 


\icdpf ^(jjcd'j 


■ p cd \ p ] vw 


dx 


F 


RC + 


Ki 

1 


Kb 


” Cd Pe{ ( t >Cd ) 


-p Cd IP 


vwdx — 


d x p 


RC + 


” cd p' e (<p cd ) 


- P cd lP 


vwdx 


Ke 


dr.V 


cd 


RC 


„cd 


i p' e {<j) cd ) 


vwdx + 


M 


/c 7 

RC + 


" cd p' e (4> cd ) 


■ p cd I p I vw 


dx 


ic 8 


K 9 


pd x <t>' 


cd 


VV { 


cd 


RC H - 


l 


\W cd p' e (4> cd ) 


-p cd ,p 


vwdx. 


/Cio 


(4.8) 


We now turn to compute Ki (1 < l < 10) term by term. For the delicate term 1C i, it can be rewritten as 

1 


JCi = — f ipdt (RC ~ P cd P) vwdx — f ipdt ( 
J R + J R + V 

= — (7 — 1 ) / ipvwdt 

J Rj. 


R 

7 - 1 ' 


*+ 

„cd, 


" cd p' e ((p cd ) 


p ) vwdx 


C + p c p ) dx + 7 / 1 pvwp d x ipdx 

S.+ 


+1 L \ L dd * 2) (^WFy) dx - L * vSt (-^hn) 


v cd p' e (cp cd ) 

= ( 7 - 1 ) f ipw (RC -p cd p) (d x u cd + d x ip) dx+n(pi - 1 ) f 
J Rj. J Rj 


vv cd p' e (4> cd ) 

v cd d x C - d x e cd p 


vwdx 


,,cd 


d x (ipvw) dx 


K 1,1 


K 1,2 


(7 — 1 ) j ipvwGdx — C f p cd d x vwip 2 dx — C f d x p cd vwip 2 dx + 7 f ipvwd t p cd pdx 
J r + 2 7k + 2 J R+ J R+ 


(4.9) 


K 1 ,: 


Ki,. 


Ki,i 


Ki, ( 


1 - V v ) ^ (i 33 ^) 


/Ci,7 


tppd t 


vv cd p' e (<p cd ) 


vwdx , 


/Cl, 8 
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where in the third identity we have used 


R 


7-1 


dtC + P cd dtV = - RC pC V ( d x u cd + O x ,/>) + nd x 


v cd d x ( - d x O cd <p 


n cd 


+ G, 


which is derived from m and m ■ 

Since p' e {(j) cd ) < 0 according to the assumption (.4) 2 , (14.91) further implies 

0 < -1 ' i”"”) * 2 ^ 2dx = + 1>»- < 4 - 10 > 

To compute /Cyz (1 < l < 8) and K-i (2 < l < 10), by applying (12.191) . (14.41) . Cauchy-Schwarz’s inequality, 
Sobolev’s inequality (12.201) . Young’s inequality and Lemmas 14.21 and fOl we directly address the following 
estimates: 

|ACi,i| + |/Ci, 6 | + |£i, 8 | + \IC 2 \ <C f \i>w\(\d x u cd \ + |d t 0 cd |)(|C| + M )dx + C f |#>0*V>|(|CI + M )dx 

JR + JR + 

<C\\wd t 0 cd \\ L ~ || II 2 + CS\\d x ^W 2 + 5 [ w 2 i/) 2 [ip, C,] 2 dx 

0 Jr + 

<C6\\d x ^\\ 2 + C6e 0 (1 + t)~i + ^\\w\\U IH| 2 |lb, C]f 
<C6\\d x ^\\ 2 + CSe 0 (l + t)-i, 

|/Ci )2 | + |ACs| <C f (\d x <p\ + |9 X CI + \d x 9 cd ip |) {\d x ipw\ + \ip(d x 0 cd ) 2 \ + \ipwd x v\) dx 

JR+ 

<C6\\d x [<p,M\\ 2 + C8(l+t)-i, 

|/Ci i3 | <C f \ipw ((d x u cd ) 2 + (d x ip) 2 ) I dx < C f i/j 2 w 2 (d x u cd ) 2 dx + C f (d x u cd ) 2 dx + C5\\d x i/.’\\ 2 
JR+ 

<CS\\d x ip\\ 2 + C5(l + t)-i, 

|/Ci, 4 | + |/Ci ; 5 | + |/Ci, 7 1 < C [ \d x 9 cd wip 2 \dx + C f \d x (pwi/j 2 \dx 


<CS [ {d x 9 cd ) 2 ^ 2 dx + C6\\d x [^,p}\\ 2 + CS(l + t)~ 2 , 

J R+ 

\ic 3 \ < c [ IV’dCI + W\)dtw\ dx < C5e 0 {l + t)~i, 

J M_)_ 

I/C 4 I + |/Cio| < C f | d x 9 cd w (C 2 + <p 2 ) | dx + C f | d x ipw (( 2 + ip 2 ) \ dx 
J R + J R + 

<C5(1 + t)~ 1/2 [ \d x 9 cd \\x\(\\d x C\\ 2 + \\d x y\\ 2 + \y(Q,t)\ 2 )dx 

JR- 1_ 

+ C6\\d x pX\\ 2 + j\\w\\l~ (IICIIIM + IMIIM) IIMII 2 
<C5\\po\\ 2 H ie <* * + C5\\d x [tp, C ]|| 2 + CS{1 + t) 2 , 

|/c 6 | <C6 2 (l + t)-i\\F\\ Ll (IICII00 + Halloo) < C6 2 (l + t)-i (||C||^R# + IMI*IIM*) 
<C8\\d x [p, C]|| 2 + C5(l + 1) 3 , 

|/c 7 | <c f la^o-dd + MH dx < C 5 (i + t)~* (IICII00 + Halloo) R^IIIMl 

J M_|_ 

<cs(i + t)-i (iKII^RCII^ + IMIMI 

<C8\\d x [y, C]|| 2 + C8(l +1)~ 2 , 
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|ACg| < C / \d x O cd w\\a\{\C,\ + \y\)dx <C5\\ip 0 \\ 2 m e » + C8\\d x [<p, <, a 


|/Cg| <C / [\dl<j\ + \d x ad x v\ + \a 2 \ + \d x v cd d x v\] [flc^CI + \d x <p\)\w\ + (|C| + \p\){{d x 0 cd ) 2 + \d x vw\)] dx 

J ]R_|_ 

<C8\\d x [(pJ,a,d x a)\\ 2 + C8(l + t)~i. 

We now plug the above estimates for JCij with 1 < l < 8 into (14.1011 to obtain 


/Ci - 


- ip 2 (d x d cd ) 2 dx 


2 v cd p' e (ft cd ) 2 

<^11^^,^,C]|| 2 + C<5(l+t)-i+CT [ (, d x 6 cd ) 2 il> 2 dx. 

JR * 


(4.11) 


Next by substituting the estimates for /C; (2 < l < 10) and (14.1111 into (14.811 and integrating the resulting 
equality with respect to time over [0,T], one has 


- / 'ft 

J M._(_ 


RC- 


i 


cd 


VV™ p' e {(f> cd ) 

RC, + 


vwdx 


1 


VV cd p' e (4> cd ) 
T 


-P cd \v 


1 2 


2 v cd p' e ((j) cd ) 2 

<C8 + C8e 0 + C5\\v 0 f H1 +C8 f \\d x [p, ft, C, d x a\\ 2 dt, 

Jo 


+ ^-p cd v - CS ) ip 2 } (, d x e cd fdxdt 


for suitably small <5 > 0 and e 0 > 0. 
Let us now define 


£M,C) = - [ 

J R + L \ vv Pe 


{ft 


,cd\ -P^IV 


vwdx , 


then (14.31) follows from (14.121) and (14.131) , this ends the proof of Lemma 14.31 


(4.12) 


(4.13) 


□ 


Finally we give the detailed proof of (12.281) . 

Proof of (12.281) . Taking the inner product of d* \2A\ with d t a with respect to x over R + , one has 


..cd 


..cd 


d t { — dzad t adx + / — d t dzad t crdx - / d t —rrd x v d x ad t adx- / -^-d x vd t d x ad t adx 


.cd 


..cd 


— ( d t pd t adx + f d t v (l — v cd p e {a + ft cd )) d t crdx 
+ f vd t (l - v cd p e {a + (j) cd )) dtcrdx - f d t (v cd d ; ^ x ^ 

Jr+ Jr+ V 


3a 


dtcrdx. 


(4.14) 


We turn our attention first to Ji (1 < l < 4) which can not be directly controlled. Since a(0,t) = 
cr(+oo,t) = 0, by integration by parts and using the cancellation, we find 

r v cd f F) v cd 

J -2 + Ji = - / - d t d x crd t d x adx - / —- d t d x ad t crdx, (4.15) 

J R+ V J R+ V 

and 

r f). v cd „ r v cd r B*v cd 

Ji+Jz= / - d 2 crd t .crdx + / — T d t vd x ad x d t adx - / d x vd x ad t adx. (4.16) 

J r + v J R+ v J R+ v 
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On the other hand, similar to (12.111) . one has 


1 -v cd p e {a + <j> cd ) 


= -v cd p' e {(j) cd )(j-v cd 



~ Q)dQ, 


Jo 


(4.17) 


and moreover 

d t Jo ~ d t cra + d t v cd a + d t v cd a 2 . (4-18) 

Substituting (14.151) . (14.161) . (14.171) and (14.181) into (14.141) and applying (12.11) and (11.71) . we deduce 

,.cd 


[ — \d t d x a\ 2 dx- [ vv cd p' e ((f> cd ) \dt.a\ 2 dx 
JKj. v J Rj. 


< 

f r) v cd 

/ —- d t d x ad t crdx 

+ 

f ^ alaa,adx 

+ 

/' v cd 

/ d t vd x crd x d t adx 


J R + V 


J R+ V 


J R + V 2 


d t v‘ 


-dxvdxcrdtcrdx 


+ C f \d x i/jdtcr\dx + C f \d x ucrdtcr\dx 
*/R+ 

+ C ( \dtv cd adtcr\dx + C f \dtaadtcr\dx + C f \d t v cd a 2 d t cr\dx 
*/R_|_ J R_|_ «/R_|_ 

+ C [ \d t d 2 v cd d t a\dx + C [ \d t v cd d 2 x v cd d t a\dx + C [ \d t v cd d x v cd d x vd t a\dx 
t/R-l- R_|_ */R_|_ 

+ C f \d x v cd d 2 ud t (j\dx, 

J R + 

which yields (12.28L according to Cauchy-Schwarz’s inequality, (12.191) and Lemma [4.11 This completes 
the proof of (12.281) . 
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